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The role of the dilaton field and its coupling to matter may result to a dilution of Dark Matter 
(DM ) relic densities. This is to be contrasted with quintessence scenarios in which relic densities are 
augmented, due to modification of the expansion rate, since Universe is not radiation dominated at 
DM decoupling. Dilaton field, besides this, affects relic densities through its coupling to dust which 
tends to decrease relic abundances. Thus two separate mechanisms compete each other resulting, 
in general, to a decrease of the relic density. This feature may be welcome and can rescue the 
situation if Direct Dark Matter experiments point towards small neutralino-nucleon cross sections, 
implying small neutralino annihilation rates and hence large relic densities, at least in the popular 
supersymmetric scenarios. In the presence of a diluting mechanism both experimental constraints 
can be met. The role of the dilaton for this mechanism has been studied in the context of the 
non-critical string theory but in this work we follow a rather general approach assuming that the 
dilaton dominates only at early eras long before Big Bang Nucleosynthesis. 
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I. INTRODUCTION 

The nature and origin of the Dark matter (DM) of our Universe is one of the big mysteries of modern 
Cosmology whose resolution is still pending. Analyzing the data cumulated from various observations 
over the past twelve years it is found that 96% of the Universe energy budget today consists of unknown 
entities, 23% of which is Dark Matter and 73% Dark Energy (DE), or vacuum energy, which is responsible 
for the current acceleration of the Universe. These data include, among other, observations of the 
Universe acceleration, using type-la supernovae [l| , measurements of cosmic microwave background 0, Q 
anisotropics, baryon oscillation Q and weak lensing data [H. The aforementioned results follow from 
best-fit analyses of various astrophysical data to the Standard Cosmological Model (ACDM) which can 
successfully describe the evolution of our Universe. The model is based on a Friedmann-Robertson- 
Walker (FRW) cosmology, involving cold DM, at a percentage 23% , baryonic matter at 4 % and a 
positive cosmological constant A > which is put in an ad-hoc manner in an attempt to describe the 
vacuum energy density. 

Supersymmetry provides one of the leading DM candidates, the neutralino which is still lacking ex- 
perimental verification. Its thermal abundance, calculated in the context of the simplest supersymmetry 
models (minimal supersymmetric model embedded in minimal supergravity is severely restricted by 
cosmic microwave background data. In the near future incorporating data from collider experiments, 
such as the LHC Q, these models may be possibly ruled out. However the existence of scalar fields in 
the primordial Universe, that contribute to the energy density, may play a dramatic role and upset the 
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whole scenery. The quintessence @ has been invoked in an attempt to explain the vacuum energy, in 
the sense that the energy it carries today is the vacuum energy measured in astrophysical observations. 
Its existence affects the relic abundances if Universe is not radiation dominated during DM decoupling. 
In fact DM relic density is predicted to be enhanced @, and in some cases this enhancements reaches 
~ 10^ or so 'T^, for a review see for instance In general, modifications of the expansion rate and 
departures from the standard cosmological scenarios may have dramatic consequences for the DM relic 
density [l^, [3 and the observed amount of DM puts constraints on possible modifications of the Uni- 
verse expansion at early eras A particular class is the tracking quintessence scenario in which the 
quintessence field is in a kination-dominated phase at early eras [Ts ]. In this context the predictions for 
the gravitino and axino DM are considered in [l^ while in the predictions for the neutralino DM 
relic, in the popular supersymmetric schemes, is discussed in the light of the constraints arising from the 
observed e^-spectrum by PAMELA [l^ and Fermi-LAT observations jT9| . 

In some string inspired scenarios, with a time-dependent dilaton-0 sources po| . whose evolution is 
dictated by non-equilibrium string dynamics pT( . the amount of thermal neutralino relic abundance is 
diluted by factor of 0(10), relative to that calculated within the ACDM-minimal supergravity cosmology 
and such models are found to survive the stringent tests of LHC ^251 . The dilution is due to the appearance 
of a friction-like term, on the right-hand side of the ap pro priate Boltzmann equation. This term plays 
also a significant role in other considerations studied in |23| . 

In this paper we argue that the mechanism for the dilution of DM relic abundances is more general 
and can hold in other instances too, having its basis on more general features of the dilaton dynamics 
prevailing in early eras, independently of the non-criticality of the underlying string theory. In order to 
model the dilaton behavior we assume the existence of exponential-type potentials occurring in a wide 
class of quintessence scenarios, supergravity models, or arise in string theories from quantum corrections. 
The presence of such a dilaton field, that dominates over radiation long before nucleosynthesis, affects the 
predictions for relic abundances in a dramatic way. In fact the conventional calculations get smaller by 
factors as small as 0(10^^) , in some cases, allowing therefore for smaller annihilation cross sections in 
the popular supersymmetric schemes employed in literature. This may alter the potential of discovering 
supersymmertry at collider experiments since the parameter space allowed by the cosmological data moves 
to regions that would be otherwise forbidden. As far as other ways of discovering DM are concerned ( 
for reviews see the small cross sections required to explaine the cosmological data may affect the 

predictions for direct [1^ and indirect il8„ J^^J DM searches ( for a review see ). 



II. SETTING UP THE MODEL 



Omitting radiation and matter contributions the equations of motion for a time-dependent dilaton are 

4) + 3H(j) + v'{(t)) = 
= ^ + U(0) 

2ij = -ig^+P4,) = -0' (1) 

In these equations the field (p is dimensionless and the potential carries dimension mass^. The first of 
these equations is not independent but is derived from the other two. In order to model the dependence 
of the dilaton as a function of In a, where a(t) is the cosmic scale factor, we assume a linear in In a 
form during early eras, which can follow from exponential-type potentials V ^ e~'"^ . Such potentials 
are inspired by quintessence scenarios and they can also occur in string theories as perturbative or non- 
perturbative corrections. The dilaton is then given by 

^ = cln (^-^^ + (j>j , (2) 

where c is a constant and a/ = a(tj) is the cosmic scale factor at the maximal reheating temperature 
reached after inflation, denoted by Tj , which occurred at the time tj . This holds in epochs t < tx in- 
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which dilaton dominates over radiation and matter, during DM decoupHng which occured earher than 
BBN. We know that BBN took place when In {ubbn /O'o) — —22.5 , corresponding to Tbbn — 1 MeV , 
and DM decouphng occured at a temperature between Tjj m — 5 — 20 GeV , as dictated by interpreting 
DM to have supersymmetric nature, corresponding to a value of ln{aDM/ao) between ~ —31.5 and 
~ —33.0 . At times tx radiation also starts contributing to the energy-matter density and at BBN must 
overwhelm dilaton's energy. Therefore a reasonable region for which ^ holds is set by a < ax with 
ln{ax/cio) ~ — 25 or smaller. 

Beyond tx the dilaton is assumed to receive an almost constant value. The constancy of (f) when hadrons 
are non-relativistic is rather mandatory if we do not want the diluting mechanism to affect the abundances 
of the known hadrons and especially nucleons. This pushes the bound on ax , defined earlier, to even 
lower values. In fact the couplings of dilaton to matter density is through the appearance of dissipative 
terms ~ {gm — 3pm) (j^ i which modifies the continuity equation for matter, and such terms are vanishing 
when hadrons are relativistic, that is at temperatures higher than about ~ 1 GeV, coresponding to 
ln{ah/ao) ^ — 30 . Below Th however hadrons are non-relativistic and dilaton couples to hadrons 
as ^ Qm 4' ■ Therefore in this temperature regime the dilaton has to be almost constant in order to 
suppress its coupling to hadronic matter. A reasonable value is = ^qcd , with Aqcd — 260 MeV 
the characteristic QCD scale, which pushes the bound set on ax to ln{ax/ao) ~ —28.4 , although 
larger values for Th , corresponding to smaller values of In {ax/ao) , are not excluded. Such values for 
Tfi are within the range that the coupling of the dilaton to supersymmetric matter is non- vanishing and 
this may have dramatic effects for the DM relic abundances as we shall see ^ . 

For t < tx the time derivative of (f> is related to the expansion rate by <j> = cH and when this is 
plugged into the third of equations ([T]) it can be solved for the expansion rate H yielding 

= H-' + '^it-tj) . (3) 

In this and the following equations the subscript I denotes quantities evaluated at tj . For the expansion 
rate solving H = a/a , and using (jS]), we get 



fl/^r^(^-^/) + l (4) 



and therefore the time-dependence of the dilaton in this era is 

4^ = hn(^^^{t-ti) + l^+^i . (5) 

Knowing the dilaton and the Hubble rate from the second of Eqs. ([1]) the form of the potential can be 
derived 

^(^) = i^-^) Hfe-^^^-^') . (6) 



This holds in the region where the dilaton depends linearly on Ina , as in shown in Eq. ([2]), that 
is for values of the cosmic scale factor a < ax ■ If the maximum reheating temperature attained is 
Tj ~ 10^ GeV the value of the cosmic scale factor at tj is ln{ai/ao) ~ —50.86 where ao — a{ttoday) 
denotes its value today ^ . Therefore the region of applicability of Eq. ^ is for values of a satisfying 



We are aware of the fact that a constant dilaton in this range cannot account for a change Aa/a ~ 10~^ over 
cosmological time scales of the fine structure constant. Interpreting the constancy of the dilaton as small quantum 
fluctuations At/) << 1 it puts a lower limit on the couplings of dilaton to matter approaching the capability of Eotvos- 
like experiments. 

We have in mind a supersymmetric model with the MSSM content whose sparticle mass spectrum is in the TeV range. 
Under these circumstances at the temperature Tj supersymmetric as well as SM particles are all relativistic and the 
effective number of degrees of freedom is Se// = 228.75, independently of the precise sparticle mass spectrum. This 
entails to the value In (aj /ag) ~ —50.86 quoted above. 
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In (a/aj) < B where B is set hy B ^ ln{ax /ao) + 50.86 , that is a number ^ 20 or so depending on the 
value of ax- At the end point ax, f^or which In (ax/aj) = B , the potential is exponentially suppressed 

V ^exp[-c{(j)x~M]-expi~Bc^) (7) 

independently of the sign of the constant c, provided the value of |c| is not exceedingly small. Since we have 
in mind a positive potential which drops as Universe expands the constant c should be bounded by < 6 
as is evident from (|6]) . As we shall see shortly, dominance of the dilaton energy over radiation is achieved 
for > 4 and therefore the value of the potential at ax, given by Eq. ([T]), is very much suppressed 
before Nucleosynthesis, that is for times earlier than Ibbn corresponding to ln{aBBN/o-o) ——22.5. 

The ratio of the kinetic to the potential energy of the dilaton field in the regime t < tx is constant. 
This follows from the second of Eq. ([1]) and the fact that (p — cH , due to Eq. In fact 

=4-1 (8) 




This yields a ratio of dilaton to radiation energy density given by 

VW ■ (9) 



2 1 / \ 4 

nip 1 (a 



Pr pS? 1 - c2/6 \ao 

In this we have reinstated dimensions and hatted densities carry dimension energy'*. The zero subscripts 
denote the corresponding quantities today. Eq. (|9]) can be also cast in the form 

9 

p^ _ 3H^ml ( Hf\ f ai\ , „ , ^^^^ 



Pr P°r \HqJ \ao J \ai 

if one expresses the potential ^ in terms of the cosmic scale factor a . Since a > aj this ratio decreases 
for values of 4 < . 

The behavior of the dilaton and the potential as functions of ln{ a/ao ) are shown in Fig. [T] for a 
particular values of the slope c, in the range 4 < < 6, and ln{ax/ao) = —28.4. Without loss of 
generality the value (px of dilaton at the end of the dilaton-dominance period has been taken vanishing. 
Actually physics results depend on the difference (/> — 0/ so a non- vanishing value for (/)x corresponds to 
a different initial condition </)/ for the dilaton field. The ratio of the potential energy, at the end of the 
dilaton-dominated era, to the same energy at reheating temperature drops by at least fourty orders of 
magnitude. 

Since dilaton energy dominates over radiation energy in this regime the value of the Hubble rate at 
reheating temperature is constrained by Eq. (jlOp . In order to quantify this suppose the dilaton to 
radiation energy density, at a given reheating temperature Tj , is 



Pr 



IQP (11) 



then Eq. ( [TO]) yields 



Hi/H„ = X 1042+P/2 (12) 



if Tj — 10^ GeV . If inflation is responsible for the generation of the power spectrum of the curvature 
scalar Ps and tensor Pt perturbations then an upper bound on the inflationary potential, and hence on 
the corresponding Hubble rate at the end of inflation Hi , can be derived 29] . Assuming that tensor per- 
turbations are small in comparison with the scalar ones the bound imposed on Hj is Hj < —j= mp Ps' 

which in turn yields Hj < 2.65 x lO^'^GeV^. This entails to the following upper bound for the ratio 
Hi /Ho, 

1 24 

Hi /Ho < — X 10^^ . 
ho 
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FIG. 1: The dilaton (left) and the potential (right) as functions of ln{a/ao) . The ratio of dilaton's kinetic 
energy to its potential energy density stay constant and thus the fast exponential drop-off of the dilaton potential 
indicates that the dilaton to radiation energy ratio has been suppressed before Nucleosynthesis. 



Then on account of (|12l) an upper bound on p is derived, 

p < 28.5 (13) 

This merely indicates that the ratio of the dilaton to radiation energy density (1111) can be indeed large 
for reasonable values of the initial conditions set at Tj , consistent with the bounds put on Hj. 
From Eq. (jlOp the ratio Pcj,/Pr can be expressed in terms of its value at Tj through 

P± = ( P±] 

Pr \Pr) 

and if this ratio at T/ is as given in Eq. (|lll) , the corresponding ratio at the end of the dilaton dominance 
period is 

= IQP feV ' = lQV+N{i-c') . (15) 
X J 

In this equation the constant N is given by 

N = log{e){b-r) =OAMA{b-r) (16) 

with h and r given by 6 = ln{ax/cLo) and r = ln{ai/ao) respectively. The first is an input while 
the second is determined by the input value of the reheating temperature Ti . For instance for T/ = 
(T 10^ GeV the constant r is given by r = —50.86 — Ina in a model with the content of the MSSM and 
mass spectrum much lighter than T/. is a number between 8.20 and 11.16, when ln[ax / clq) is 

taken within the range —32.0 25.0 , and the value of a specifying the reheating temperature is of the 

order or unity or so. At the end of the dilaton-dominance era the dilaton to radiation energy ratio drops 
from its initial value to 

= 10"' . (17) 

X 
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The power p' must be smaller than p, but still large enough to guarantee that the above ratio is much 
larger than unity, so that the bulk of the total energy is carried by the dilaton in the regime aj < a < ax ■ 
From Eqs. p3|) and (IT71) we deduce that 

c^=4 + ^ (18) 

and therefore the drop-off of the ratio p pr yields that is larger than 4. Combined with its upper 
bound < 6 discussed earlier one concludes that lies in the rather narrow range 4 < < 6. One 
can also utilize the relation (|18p . in combination with the bound < 6 , to derive the following bound 
on p — p' 

p-p' < 2N . 

Since A'^ is a number of order ~ 10 the above upper bound set on p — p' leaves much room for values 
of the powers p, p' to guarantee that the dilaton energy indeed overwhelms radiation energy in the whole 
regime aj < a < ax , as is assumed in this scenario. 

III. DILUTION OF DM ABUNDANCES 

Concerning the calculation of relic densities, omitting the collision terms, the energy-matter density 
obeys the following equation 

^+3H{p + p)-±{p-3p) = (19) 

where the last term is the coupling of the dilaton to the density ^ . 

Obviously for radiation the last term drops and hence when matter is relativistic this term is absent. 
This holds at a temperature T » m where m is the mass of the particle under consideration. Including 
the collision terms into Eq. (1191) entails to the following equation for the number density 



^ + 3Hn+ <va> {n^ - ) _ A „ ^ Q . (20) 
dt ' \/2 

This is suitable for describing the evolution of the number density during periods for which the particle 
is non-relativistic and pressure practically vanishes. During eras in which the particle is relativistic the 
last term drops, since the last term in Eq. (IT51) drops too, and in this case (UHl) receives the well-known 
form of Boltzmann equation, [23 | . 

For the number to entropy density ratio, Y = n/s , Eq. ((20|) takes on the form 



— 1/2 

— = ax) m (va) ( -^9eff ) (^^ + 3 ^) " K) + ^i^) Y . (21) 



In this X stands for x = T/m where T is the photon gas temperature related to the radiation density 
Pr , which includes all relativistic particles at a given epoch, through 



^ The division by y/2 is due to the normahzation of the dilaton whose kinetic energy appears as p^*" = ^ , see the second 
ofEq. (dJ. 
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FIG. 2: The LSP Dark Matter number density to entropy density ratio q — 



as function of a; = 



particular supergavity model. The values of ^, 5* denote the status of the ^-factor and the source respectively ( / 
for open, x for switched-off ). For comparison the corresponding equilibrium density qo{x) has been also drawn. 



and Gat is Newton's constant. In ([2T|) the quantity h stands for the entropic degrees of freedom related 
to the entropy density through s — 2tt^ h[T) / 45 . 
The prefactor S_(x) appearing in Eq (PT|) is given by 



Pr 



-1/2 



while the source S{x) , in the same equation, is 

S{x) = 



^ X dh\ 
3h dx J 



(22) 



(23) 



In the expression for ^(x) above, Pr,Pm and are the radiation, matter and dilaton energy-densities 
respectively. Recall that we use a dilaton density having dimensions , see Eq. [1] Note that no 
cosmological-term contributes to Eq. ([22l) since such a term is absent at DM decoupling and long after it. 
In conventional treatments the prefactor ^(x) is unity since the DM freeze-out is assumed to take place 
in the radiation dominated era. However in the presence of the dilaton energy term this is smaller than 
unity and the density Y decreases slower, as temperature drops, than in the conventional cases where 
£,{x) — 1 . In the source term S{x) , given by Eq. (1^51) . the quantity (/>' is the derivative of the dilaton 
with respect ln{a/aa) and if this is negative then the source acts in the opposite direction than S,{x) 
tending to decrease the density Y faster as x decreases. 

It should be noted that for simplicity we have assumed the lowest order, in a', contributions to the 
form factors e~^^'^^ and Z{(j)) associated with the scalar curvature R and dilaton kinetic terms of 
the effective action, in the string frame, and hence the simple expressions for the ^-dependent terms of 
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q{x) = n/T^/h 

le - 005 



le - 006 



le - 007 
le - 008 
le - 009 



le - 010 



le - Oil 
le - 012 



le - 013 

le - 014 

0.01 0.0001 le - 006 le - 008 le - 010 le - 012 le - 014 
X = T/m 

FIG. 3: The same as in figure [2] with x — — from 0.1 to values corresponding to CMB temperature today. 



Eqs (fT9|) and (pO| . Also the dilatonic charge has been assumed vanishing. However the couplings of 
the dilaton to matter may evolve in time with the dilaton itself and depend on the particle species in 
a non-universal way. Therefore other options are available which in the string theory arise from loop 
corrections or non-perturbative string effects [10, IMj . In such cases the coupling of matter to cj) in the 
continuity equation Eq. (|19p. which is mainly controlled by ip'iff') = dip/dcj), is not a constant. Besides 
there is an additional contribution to the continuity equation which depends on the dilatonic charge, if 
the latter is assumed non- vanishing [sol - fs^ . Including these effects will give rise to modified (/(-dependent 
terms in Eqs. (fT9)) and (l20l) resulting to a source S{x) in Eq. ([23)1 which is multiplied by ^ '0'(0), 
provided the dilatonic charge of Dark Matter is taken vanishing. This will still tend to decrease the 
density Y, as x decreases, if tp'{4>) > in the regime following Dark Matter decoupling. In particular 
if '0'(0) > 1 the dilution caused to the relic density is enhanced, in comparison with that caused by the 
source term as it appears in Eq. ([23]) . or gets smaller if ?/''('/') < 1- Certainly in order to further study 
the effects of this term one needs a better understanding of how to handle the corrections to ip{(l)) arising 
from the underlying string dynamics. For definiteness in this work, and in order to quantify the effect of 
the dilution of the abundance of Dark Matter, we assume a gravi-dilaton effective action in the lowest 
order in the string slope a'. 

The effects of the presence of the factors £,{x) and S{x), as given by Eqs (1^ ) and ([Ml), is shown 
in figures [5] and [3] where we plot the density as function of the temperature, actually x — T/ttilsPj for 
particular SUSY inputs. Eq ([?!]) has been integrated numerically which yields more accurate results than 
the approximate solutions employed in (20| . The displayed figures correspond to a supergravity model 
with inputs given by mo = 1100.0 GeV, M1/2 = 1200.0 GeV and Aq = GeV . We have taken tanP = 40 
and the parameter /i is taken positive, /i > 0. The value of 6 = ln{ax/ao) , setting the onset of the epoch 
after which the dilaton is constant (0 = 0), has been taken —28.4 corresponding to a temperature 
Aqcd = 260 Gey as we have already discussed. For the particular SUSY inputs the LSP Bino has a 
mass niLSP — 527.2 GeV and the point b = —28.4 corresponds on the x axis to a value x ~ 0.0005. For 
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comparison, except the ordinary case scenario, where both the source S and the ^— factor are absent ( 
red soHd hne ), the cases where both ^ and S are open ( green dashed-dotted Une ), or when only the ^— 
factor is present ( blue short-dashed line ) are also shown. The very thin dashed line, that rapidly drops, 
is the corresponding equilibrium density. In the most interesting case, that both terms are switched on, 
the density is monotonically decreasing after decoupling due to the appearance of the source term. The 
rapid change around x ~ 0.0005, corresponding to 6 = —28.4 where dilaton reaches its constancy, are 
shown in figures [5] and [31 In the specific example shown the relic density is diluted by a factor of ~ 50, as 
can be seen by comparing today's density values for the conventional case ( red solid line ) and the case 
where both ^— factor and the source term are present ( green dashed line ). In the first case the relic 
density predicted is ^lsp = 6.059 while in the second case the relic density is considerably reduced 
faUing in the WMAP allowed range riiSp/iQ — 0.1116. In general, for given b = ln{ax/ao) , one can 
obtain reduction factors in the range 0(5 — 50) , the smaller (larger) corresponding to lighter (heavier) 
neutralino masses. 



IV. CONCLUSIONS 



In this letter we have shown that the dilaton dynamics during early eras, long before Nucleosynthesis, 
in conjunction with its coupling to Dark Matter may have dramatic consequences for the predicted Dark 
Matter relic density. Modeling the dilaton evolution to be that dictated by exponential type potentials 
V ~ e~^'^ , occuring in quintessence scenarios and string theory, the ordinary predicted DM density 
may be diluted by large factors ranging from 0(5) to 0(50) . This dilution mechanism is consistent 
with the absence of dilaton couplings to ordinary matter ( hadrons ), in the continuity equations, but it 
affects DM relics since dilaton dominates over radiation during and after DM decoupling. This allows for 
LSP annihilation cross sections, in the popular supersymmetric schemes, that are smaller by an order of 
magnitude or more. This however may imply smaller inelastic cross sections of the neutralino LSP with 
nucleons put ting farther the potential of discovering supersymmetric DM at proposed direct detection 
experiments [26|. As far as indirect detection experiments are concerned, indirect searches of Dark Matter 
through antimatter production has stirred much interest the last three years. The PAMELA data [l^ 
in combination with those provided by Fermi-LAT fTo*! and HESS ^Tf\ may be conditionally explained 
as DM annihilations in the galactic halo that generates the produced antiparticle flux [33|]. In the case 
under consideration the smaller annihilation cross sections, required to satisfy WMAP data, makes even 
harder the possibility that antimatter fluxes observed in the cosmic ray are relevant to annihilation of the 
neutralino LSP in the galactic halo. In conventional supersymmetric models the situation can be rescued 
at the cost of considering large boost factors. However even in this case DM annihilation looks a rather 
remote explanation for interpreting the aforementioned data and other more conservative explanations 
exist, as for instance antimatter produced by pulsars or supernovae etc [3^ . 

A complete phenomenological study of supersymmetric models addressing all these issues is in progress 
and the results will appear in a forthcoming publication (35| . 
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